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ABSTRACT 
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' In light of recent observations of spectra of the quasars SDSS 1030+0524 (Becker 

et al. 2001) and SDSS 1044-0125 (Djorgovski et al. 2001), we study the observational 
■ signatures of different stages of the reionization epoch. During the initial, pre-overlap 

stage, the hydrogen throughout the universe is neutral except for isolated H II bubbles 
ff^ \ due to individual ionizing sources. The central stage of reionization is the overlap stage, 

\ during which the individual H II bubbles overlap each other and reionize the low-density 

' gas which takes up most of the volume of the universe. Some neutral hydrogen remains 

in dense clumps which are then slowly reionized during post-overlap. We show that 
both of the recent observations are consistent with the post-overlap stage. Becker et 
al. (2001) may have observed the universe in the era before the end of overlap, but a 
Q . conclusive proof of this requires the observation of similar intervals of Gunn- Peterson 



absorption at the same redshift along several additional lines of sight. 
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1. Introduction 



Until recently, the highest redshift detected for a quasar was z = 5.8 (Fan et al. 2000). Although 
the spectrum showed strong absorption short-ward of the Lyman-a line of the quasar, the detection 
of transmitted flux in this spectrum was taken to imply that reionization was complete by z ~ 6, 
since even a small neutral fraction in the IGM would have produced total absorption, a fact first 
pointed out by Gunn & Peterson (1965). The most natural explanation for reionization is photo- 
ionizing radiation produced by an early generation of stars and quasars; recent calculations of 
structure formation in cold dark matter (CDM) models find that reionization should naturally 
occur at z ~ 6-12 (Haiman & Loeb 1997, 1998; Gnedin & Ostriker 1997; Chiu & Ostriker 2000; 
Gnedin 2000, and for comprehensive reviews on the subject of reionization and the first galaxies, 
see Barkana & Loeb 2001; Loeb &; Barkana 2001). 
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The reionization of hydrogen is expected to involve several distinct stages. The initial, "pre- 
overlap" stage (using the terminology of Gnedin 2000) consists of individual ionizing sources turning 
on and ionizing their surroundings. The radiation from the first galaxies must make its way through 
the surrounding gas inside the host halo, then through the high-density region which typically 
surrounds each halo. Once they emerge, the ionization fronts propagate more easily into the low- 
density voids, leaving behind pockets of neutral, high-density gas. During this period the IGM 
is a two-phase medium characterized by highly-ionized regions separated from ncTitral regions by 
sharp ionization fronts. The central, relatively rapid "overlap" phase of reionization begins when 
neighboring H II regions start to overlap. Whenever two ionized bubbles are joined, each point 
inside their common boundary becomes exposed to ionizing photons from both sources. Therefore, 
the ionizing intensity inside H II regions rises rapidly, allowing those regions to expand into high- 
density gas. This process leads to a state in which the low-density IGM is highly ionized and ionizing 
radiation reaches everywhere except for gas located inside self-shielded, high-density clouds. Some 
neutral gas does, therefore, remain in high-density structures which correspond to Lyman Limit 
systems and damped Lyct systems seen in absorption at lower rcdshifts. The ionizing intensity 
continues to grow during this "post-overlap" phase, as an increasing number of ionizing sources 
becomes visible to every point in the IGM. 

Recently, Djorgovski et al. (2001) and Becker et al. (2001) found, for the first time, broad 
regions of Gunn-Peterson absorption in quasar spectra, with high optical depth throughout each 
region. Djorgovski et al. (2001) observed the quasar discovered by Fan et al. (2000) with a higher 
spectral resolution and high signal to noise, and discovered a dark region which is ~ 5 Mpc long 
at z ~ 5.3. Becker et al. (2001) discovered the most distant source of Ught observed to date, a 
quasar at z = 6.28; the spectrum showed a dark region of ~ 15 Mpc length extending right up to 
z = 6.15, where the emission of the quasar itself may begin to have an effect. Since even a small 
neutral fraction in the IGM can produce these long stretches of strong absorption, it has thus far 
been unclear whether these observations indeed herald the empirical discovery of the overlap era; 
instead, these observations may only be approaching the end of overlap, which itself occurred at 
an even higher redshift than that probed by current observations. 

In this paper, we construct a model of the statistics of neutral and ionized regions during 
reionization, in order to interpret the recent observations. Our goal is to use the observations to 
infer the state of the IGM at z ~ 6, and specifically to determine whether the universe has indeed 
been observed during or even before the overlap stage. We focus on the observable property which 
is not found at low redshifts but is a new feature of the recent observations, namely long intervals 
of continuous strong Lja absorption. The model must be statistical, since the ionizing intensity in 
the IGM is extremely inhomogeneous during reionization and for some period afterwards. Given 
the small sample of lines of sight probed thus far, we seek to address the question of whether a 
statistical fluctuation could produce the observed Gunn-Petcrson absorption even if the entire IGM 
has already been rcionizcd at the observed redshift. Since the properties of ionizing sources at these 
high redshifts are speculative at present, having almost no observational constraints, we strive to 
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obtain model-independent conclusions. We present our model in §2 and the results in §3. We then 
confront the model results with the observations in §4. Finally, we give our conclusions in §5. 

The basic theoretical framework in which the matter content of the universe is dominated 
by CDM has recently received a major confirmation from measurements of the cosmic microwave 
background (Netterfield et al. 2001; Lcc ct al. 2001; Halverson et al. 2001). Based primarily on 
these measurements, in this paper we use cosmological parameters ^Irn = 0.35, il\ = 0.65,ili, = 0.05, 
(jg = 0.85, n = 1, and h = 0.65, where i^rn, ^A, and $7^ are the total matter, vacuum, and baryon 
densities in units of the critical density, erg is the root-mean-square amplitude of mass fluctuations 
in spheres of radius 8 Mpc, and n = 1 corresponds to a primordial scale- invariant power 
spectrum. 



2. Modeling the Statistics of Gunn-Peterson absorption 

2.1. Ionizing Sources 

In this section we describe a model which allows us to predict the observable properties of 
quasar spectra during the various stages of reionization. Roughly speaking, the question we wish 
to address is the following; if we consider spectra taken along random lines of sight through the 
universe, at a given rcdshift, then what fraction of these spectra will be similar to those that have 
actually been observed? We define the similarity of spectra in terms of the presence or absence of 
long intervals of strong Lya absorption. The absorption depends on the distribution of gas density 
in the IGM and on the ionization state of the gas, which in turn depends on the properties of 
ionizing sources. Our strategy is to construct a model which is general enough to include a wide 
range of possible scenarios, given the current paucity of direct observational constraints on, for 
example, the properties of the ionizing sources. While we study the dependence of the results on 
the various input parameters, we especially seek a model-independent probe of the reionization 
state of the IGM. In this subsection we describe our model for the properties of ionizing sources; in 
§2.2 we model the gas density and ionization state, and obtain its opacity; finally, in §2.3 we discuss 
the physical quantities that we calculate in order to compare our model with the observations. 

Galaxies and the ionizing sources within them form in halos in which gas can accumulate and 
cool. At high redshift, gas can cool efficiently in halos down to a virial temperature of ~ 10^ K or a 

circular velocity of ~ 16.5 km with atomic cooling. Cooling is possible down to V^. ~ 2 km 
with molecular hydrogen (H2) cooling, but molecular hydrogen is fragile and is expected to be photo- 
dissociated well before reionization unless mini-quasars^ make a large contribution to the ionizing 
intensity (Haiman et al. 1997, 2000). Before reionization, the IGM is cold and neutral, and these 
cooling requirements set the minimum mass for halos which can host galaxies. During reionization. 



mini-quasar is an accreting central black hole in a dwarf galaxy. 
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however, when a volume of the IGM is ionized by stars, the gas is heated to a temperature Tiqm ~ 1- 
2 X 10^ K. Wc adopt a standard temperature of Tigm = 1-5 x 10^ K, and then the hnear Jeans 
mass corresponds to a viriahzed halo with a circular velocity of 



where this value is essentially independent of redshift. Even halos well below the Jeans mass can 
pull in some gas once the dark matter collapses to the virial overdensity. For simplicity, we adopt 
a sharp cutoff associated with this suppression, at a circular velocity of Vc = Vj /2, based on the 
results of numerical simulations (Thoul & Weinberg 1996; Quinn et al. 1996; Weinberg et al. 1997; 
Navarro Sz Steinmetz 1997; Kitayama & Ikeuchi 2000). However, this pressure suppression is not 
expected to cause an immediate suppression of the cosmic star formation rate, since even after 
fresh gas infall is halted the gas already in galaxies continues to produce stars, and mergers among 
already-formed gas disks also trigger star formation. However, if the overlap phase of rcionization is 
relatively gradual then there may be time for the suppression of star formation to become significant 
even before the end of overlap. 

Once gas collects inside a halo and cools, it can collapse to high densities and form stars or 
a mini-quasar. Regardless of the nature of the source, its ability to form is determined by gas 
accretion which, in a hierarchical model of structure formation, is driven by mergers of dark matter 
halos. Therefore, in order to determine the lifetime of a typical source, we first define the age of gas 
in a given halo using the average rate of mergers which built up the halo. Based on the extended 
Press-Schechter formalism (Lacey &: Cole 1993), for a halo of mass M at redshift z, the fraction of 
the halo mass which by some higher redshift Z2 had already accumulated in halos with galaxies is 



where D{z) is the linear growth factor at redshift z, S{M) is the variance on mass scale M (defined 
using the linearly-extrapolated power spectrum at z = 0), and Minin(22) is the minimum halo mass 
for hosting a galaxy at Z2 (as determined by the cutoff Vc which was discussed above). We define 
the total age of gas in the halo as the time since redshift Z2 where Fm{z, Z2) = 0.2, so that most 
(80%) of the gas in the halo has formed stars only since then. Thus, e.g., at z = 6 the age of the 
universe is 9.1 x 10^ while the age of gas in a 3 x 10^ halo is 3.6 x 10® yr {z2 = 8.8), and the 
age of gas in a IO^^Mq halo is 5.5 x 10® yr {z2 = 12.0). Low-mass halos form out of gas which has 
recently cooled for the first time, while high-mass halos form out of gas which has already spent 
previous time inside small galaxies. We emphasize that the age wc have defined here is not the 
formation age of the halo itself, but rather it is an estimate for the total period during which the 
gas which is currently in the halo participated in star formation. However, the rate of gas infall is 
not constant, and even within the galaxy itself, the gas may not form stars at a uniform rate. The 
details involve complex astrophysics, so we account for the general possibility of bursting sources 
by adding a parameter ^, the duty cycle. We pick a standard value of C = 0.25, which implies that 




(1) 
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each source is four times brighter during a lifetime which is four times shorter than the gas infall 
age that we specified above. This value of is motivated by the fluctuating accretion rate seen in 
simulations, but in our model we allow for a wide range of possible values of 

Consider an H II region produced by a source residing in a halo of total mass M and baryon 
fraction Since recombinations within the H II bubble are relatively unimportant, we assume 

that the radius of the bubble is close to the maximum radius rmax that is simply determined by 
equating the number of ionized hydrogen atoms to the total number of ionizing photons produced 
by the source over its lifetime. We also assume for simplicity that each bubble is near its maximum 
size; this is a safe assumption for our purposes since the neutral fraction is relatively high at the 
outskirts of each H II bubble, and so the cross-section for observing a region with low optical depth 
is dominated by the central regions of the bubbles. To estimate the radius, we assume that the 
baryons in the halo are incorporated into stars with an efficiency of /star = 10%, and that the 
escape fraction for the resulting ionizing radiation is /esc = 5%, where this low escape fraction is 
understood to include the effect of recombinations in the dense gas within the source galaxy, its 
host halo, and the relatively dense immediate surroundings of the halo. If the stellar IMF is similar 
to the one measured locally (Scalo 1998) then ~ 4000 ionizing photons are produced per baryon 
in stars (for a metallicity equal to 1/20 of the solar value; Haiman, personal communication). We 
define a parameter which gives the overall number of ionizations per baryon, 

ATio^ = 4000 /star /esc • (3) 

Denoting by n^^ the present mean number density of hydrogen, the radius of the H II bubble is 

_ J_ f 3_ ^ _ / 7 \ / iVion M 0.148 \^ 

Here and throughout this paper we use proper, not comoving, distances. Note that the radius rmax 
is larger than the halo virial radius by a factor of ~ 20 that is almost independent of redshift and 
halo mass. Also note that we would obtain a similar result for the size of the H II region around 
a galaxy if we considered a mini-quasar rather than stars. This results from the high expected 
escape fraction for mini-quasars (~ 50%, e.g.. Wood Sz Loeb 2000) which, together with the high 
radiative efficiency of ~ 6%, overcomes the low efficiency (~ 0.2-0.6%, e.g., Magorrian et al. 1998) 
for incorporating the baryons in a galaxy into a central black hole. 



2.2. Opacity of the IGM 

Before the end of overlap, some large voids in the IGM arc still neutral, and they completely 
block out the flux at the Lya wavelength. However, even at a wavelength A which corresponds to 
redshifted Lya inside an H II region, the optical depth rises with redshift. There may in general be 
two sources of opacity at high redshift. First, the gas inside the H II region absorbs light strongly 
in the area of the resonance of the Lya line, and can produce a significant opacity even with a very 
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small neutral hydrogen fraction. At higher redshifts, the mean density of the universe is higher, 
the voids are less empty because structure formation is still in its infancy, and the neutral fraction 
is higher since the ionizing intensity is lower. In addition, if the H II region is itself surrounded 
by still-neutral IGM then the damping wings of the Lya line due to this neutral gas can be broad 
enough to produce a substantial opacity over the entire wavelength range corresponding to the H II 
region. Indeed, Miralda-Escude (1998) showed that the neutral IGM can block out all H II regions 
up to a diameter of ~ 1 Mpc. 

To determine the opacity, we begin with the standard (Gunn & Peterson 1965) opacity for 
absorption which includes the resonance, 

= m^cHiz) = ^-^^ ^ io:0325 J (,035 J J ' 

for neutral gas at the mean density at redshift z. Here and throughout this paper we assume the 
high-redshift form for Bubble's constant, 

H{z)^HoVn^{l + z)^ , (6) 

where Hq is Hubble's constant at z = 0. More generally, gas at relative density A = Pg/pg and 
with a neutral fraction xri produces an optical depth 

Treson = TQP a^HI A , (7) 

where the subscript stands for resonance. 

For the density distribution in the IGM within the H II bubble, we adopt the general picture 
(though not the full model) of Miralda-Escudc ct al. (2000). This picture is based on the fact that 
voids can be quickly and easily ionized while the high recombination rate of dense clumps keeps 
them neutral. Thus, a rough description of the state of the IGM is that most of the gas up to some 
fiducial overdensity is ionized, and most of the gas above that overdensity is still neutral. Although 
the remaining neutral clumps can contain a significant gas mass, they take up only a small fraction 
of the volume and do not block a significant fraction of the ionizing radiation emanating from the 
central source. Thus the neutral clumps can essentially be ignored, and the total recombination rate 
is typically not high enough to limit the size of the H II bubble. Thus we describe the IGM within 
the large H II bubble as having an average density A relative to the mean, where it is understood 
that the 'H II bubbles' are not fully ionized since they contain a small volume of dense neutral 
clumps. Since the clumps are denser than average, the gas within the bubble which is ionized is 
characterized by a A < 1. For the model of the density distribution used in Miralda-Escude et al. 
(2000), if gas at z = 6 is ionized up to a density of 10 times the cosmic mean then the ionized gas 
has a A = 0.87 and a clumping factor (i.e., mean value of density squared) of C = 1.6 times that 
for gas at the cosmic mean density. If the gas is ionized only up to a density of 2 times the cosmic 
mean then the numbers are A = 0.67 and C = 0.60. Because of these relatively low clumping 
factors, an ionizing source can easily ionize the gas up to an overdensity of ~ 10 without sacrificing 
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many photons to balance recombinations. The actual maximum overdensity of ionized gas depends 
on the detailed geometry of the overdense regions and must vary with the distance from the source. 
Indeed, the model of Miralda-Escude et al. (2000) is most accurate when the ionizing intensity in 
the IGM is uniform, and this is true only well after overlap; therefore, in this paper we instead use a 
model based on the H II bubbles of individual sources, a model which is much more accurate during 
overlap. We use the Miralda-Escude et al. (2000) picture only in the limited sense of determining 
the effective value of A. An important detail is that while the size of the region that a given source 
can ionize depends essentially only on the mean A of ionized gas, the probability of observing a 
low optical depth in the region is particularly sensitive to the presence of the lowest-density voids. 
For the density distribution of Miralda-Escude et al. (2000), only 0.9% of the mass, and 4.2% of 
the volume, is occupied by gas with A < 0.25 at z = 6. The numbers are 7.6% of the mass and 
25% of the volume, for A < 0.4. In our calculations we thus use two different effective values of 
A: wc adopt as standard values (A) = 0.8 which is used to calculate the size of each bubble, and 
At- = 0.4 which is used to calculate the optical depth of the gas. Thus, e.g., the radius of the H II 
bubble [eq. (4)] is modified to 



1 



7 \ /iVion M 0.148 0.8 Va 



r_ = 82kpc^3^J(^^^^^^J . (8) 

We assume ionization equilibrium for the hydrogen within the H II bubble. This is a safe 
assumption even though the recombination time for the IGM at z = 6 is around the Hubble time. 
Since equilibrium typically means a neutral fraction of ~ 10~^ or less then the time to reach 
equilibrium is at most 10~^ of the recombination time. The recombination rate per volume is 

A^ 

Rec rate = i'n'fj)'^ — , (9) 

where a = 1/(1 + z) and the case B recombination coefficient for hydrogen at T = 10^ K is 
aB = 2.6 X 10~^^ cm^ s~^. The ionization rate per volume is 

dN^ 1 _ _n A 
Ion rate = xhi — t"^ -; — ^ t nzr (10) 

at 47rr^ a-^ 

where a is a frequency-averaged photoionization cross section, 2 x 10^^^ cm^. We set dNj/dt = 
N^/tg, where A^'-y is the total number of ionizing photons produced (which equals the number of 
baryons in the halo times Mon) and tg is the source lifetime (see §2.1). Then the neutral fraction 
at radius r in the H II bubble is 



--^■^^^^""Irjino^ji^j (=ijP)I^I • <"> 



and the optical depth at Lya caused by this gas is therefore 
= 455 



' reson 



0.1487 VO-0211 1.3xl08yry V 20 IO^MqJ \0A 82kpc 



(12) 
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Note that the optical depth at a given radius depends on A,- (by definition of A^), but it does not 
depend on rmax (and thus on (A)), assuming that r < rmax- 

Our choice of duty cycle C affects the statistics of absorption not only through the lifetime 
ts oc if, e.g., C, = 0.25, then for each active source with a highly ionized H II bubble there are 
three 'dead' bubbles, regions ionized by sources that have already turned off before z = 6. These 
bubbles have typically been dead for ~ 2 x 10^ yr, while the recombination time (if (A) = 0.8) is 
~ 2.5 X 10^ yr. Thus, if the bubbles were initially highly ionized, at z = 6 they have a neutral 
hydrogen fraction below 10%. This fraction is high enough that the internal optical depth in 
these bubbles satisfies Treson S> 1, but the fraction is low enough that this gas does not contribute 
significant absorption through the Lya damping wings. Thus, we treat the dead bubbles as neutral 
for the purpose of internal optical depth, but as fully ionized for the purpose of damping wings and 
in the calculation of the total filling factor of ionized gas in the universe. 



2.3. Specific Predictions 

In order to analyze the importance of the observations, we look for a specific property of the 
observed spectra that can be computed from our model. Although the mean optical depth declines 
with redshift, and this quantity is straightforward to determine observationally, we prefer to focus 
on the essential feature which defines Gunn-Peterson absorption. This is not the value of the mean 
optical depth, which is already rather high even at a somewhat lower redshift; the key feature is 
instead the presence of a long stretch of absorption which is continuous, i.e., which docs not allow 
through any significant fiux at all, even when observed with a very high spectral resolution. This 
feature also captures the physical difference between the IGM before and after reionization. At 
low redshift, observed spectra are characterized by overall transmission except for a dense forest 
of relatively narrow absorption lines corresponding to the cosmic web of sheets and filaments. The 
Lya forest becomes denser with redshift, but as long as the low-density gas in the IGM is very 
highly ionized, a high-resolution spectrum should always show some transmission spikes, even when 
these spikes are rare and the mean optical depth is high. At redshifts approaching reionization, on 
the other hand, the IGM itself has a significant neutral fraction, and it produces continuous strong 
absorption over large redshift intervals. Before full reionization, parts of the IGM are neutral and 
they produce even longer intervals of total absorption. 

Prom the previous subsections, our main parameters and their default values are the redshift 
z = the minimum halo circular velocity Vc = 16.5 km/s of halos with galaxies, the relative 
density At- = 0.4 of ionized gas within the bubbles (where also (A) = 2At-), the overall source 
efficiency A'^ion = 20, and the duty cycle C, = 0.25. Por given parameter values, we calculate the 
total filling factor Qn of H II bubbles, which is simply the integral of halo number density times 
the maximum ionized volume 4mr^^-^/2> due to each halo. We emphasize that the volume taken 
up by H II bubbles includes the small volume fraction of neutral clumps which remain within each 
bubble. The end of the overlap stage of reionization corresponds to Qn = 1, because at this point 
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all the low-density gas is ionized. We allow Qu to be greater than unity in the post-overlap stage, 
where, e.g., a value of Qu = 5 means that the total number density of ionizing photons produced 
by all sources is five times greater than the mean number density of hydrogen atoms in the IGM 
(where the dense neutral clumps are not included). For direct comparison with observations, we 
compute also the mean free path Aqp for observing a region with optical depth less than a given 
f. We deal with the mean free path, a statistical quantity, since the ionizing intensity in the IGM 
is extremely inhomogeneous during reionization. At redshift z, we have 



where n{M) is the proper number density of halos with mass up to M, and erf (M) is the cross- 
section for intersecting some point with a Lya optical depth less than f along a line of sight through 
the H II bubble produced by a halo of mass M. The mass Mmin is determined by the halo with the 
minimum or, if it is bigger, by the smallest halo whose H II bubble is large enough to produce 
some point with r < f. In a stretch of length I along a random line of sight, the average intersected 
number of H II bubbles with a r < f region is Z/Aqp- Also, the chance of not intersecting any 
such bubble in the same stretch is exp[— Z/Agp]- We can similarly calculate a mean free path for 
encountering any ionized gas, regardless of the optical depth: 



where in this case Mmin is always set by the cutoff Vc, and aii{M) = Trr'^^g^^. It is also easy to show 
that a fraction Qu of a random line of sight is covered by some part of an H II bubble, so on average 
a stretch of length An intersects one bubble with a typical intersected path of length Qu^ii- 

In order to obtain Agf, we must first calculate the optical depth of every parcel of gas. To 
calculate the optical depth at each point within an H II bubble at redshift z we add Treson to the 
optical depth caused by the red damping wings of neutral IGM at lower z and the blue damping 
wings of neutral IGM at higher z. We include the blue damping wings since we are considering 
here bubbles along the line of sight to a given source (e.g., a quasar), and are not considering 
the region ionized by the quasar itself (i.e., the proximity effect). The optical depth is given by 
an integral along the line of sight of the number density of neutral hydrogen times the scattering 
cross-section of the Lya line. Consider at high redshift a stretch of IGM at the cosmic mean 
density which has a neutral fraction a^Hi and is located between redshifts Zi and Z2- On light with 
an observed wavelength A = Aq;(1 + z), this gas produces through the damping wing an optical 
depth of (Miralda-Escude 1998) 




(13) 




(14) 




(15) 



where rop is given by eq. (5) and 
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and we have assumed z < zi < Z2 oi zi < Z2 < z so that the absorption does not involve the central 
region of the Lya line, near resonance. 

The optical depth produced by an IGM with a given filling factor Qn depends on the distri- 
bution of neutral and ionized regions within the IGM. When Qu > 1 we assume that there is no 
neutral IGM and that rjamp = 0. In this case, not only is there no neutral IGM, but also most 
points in the IGM feel the radiation from more than one source. As a simple way of accounting 
approximately for the high number density of sources, for each source we consider its bubble only 
out to radius Vmax/VQu, since the more distant gas is likely to be closer to some other source, 
given that Qu > 1. Note that we are still assuming that the radiation intensity at each point is 
dominated by the nearest source; we discuss this issue further in §4. 

When Qii < 1, Tdamp depends on the exact distribution of neutral and ionized gas on both 
sides of a given H II bubble. It is not possible to analytically describe the full complexity of this 
distribution, so we use the statistical properties of this distribution in order to get a handle on 
the expected range of values of Tdamp- Specifically, we know the volume-averaged neutral fraction 
xhi = 1 — Qu, and we also know that on average a stretch of length An intersects a single bubble 
with atypical intersected path of length QuAn [see equation (14)]. We therefore adopt the following 
simple model for the IGM on each side of the H II bubble, along a given line of sight: We take 
two consecutive segments of length An each, where in each segment the central length of QiiAn is 
ionized and the rest is neutral. After these two segments, we simply take a uniform IGM with the 
neutral fraction set equal to 1 — Qu- The optical depth is only weakly dependent on the precise 
redshift at which this neutral IGM is cut off, so we fix a minimum z = A and a maximum z = 12 
(i.e., we do not consider sources observed at z > 12). By trying other possible arrangements of 
neutral gas within the IGM we find a typical systematic uncertainty of around a factor of two in 
Agp (when Qu < 1); the particular arrangement of segments that we have chosen is arbitrary, 
except that it gives results that lie in the central range of values obtained from various possible 
arrangements. 

To illustrate some numerical values, since with our standard parameters Qu = 5.0 wc consider 
here an example with z = 8, A^ion = 10 and 14 = 40 km/s, for which Qu = 0.39. In this case. 
An = 0.29 Mpc. Consider, for example, a lO^^M© halo, which produces an /"max — 0.24 Mpcj if a 
line of sight passes at a projected radius of 0.05 Mpc from the center of this bubble, the total length 
through the H II bubble is 0.46 Mpc. At the center of this line through the bubble, the optical 
depth is 11.1, made up of Treson = 5.6 and Tdamp = 2.7 and 2.8 from the red and blue damping 
wings, respectively. Note that in calculating the damping wings we have neglected the possibility 
that they could be affected by density fiuctuations or peculiar velocities in the neutral gas. The 
overall Tdamp results mainly from the integrated effect of a length of ~ 0.5 Mpc of neutral gas. At 
z = 6, the typical 1 — a density fluctuation on this scale is 0.32, so fiuctuations are likely to have 
only a moderate effect on the overall statistics, but the effect is significant and should be studied. 
The damping wings block out small bubbles and imply that typically, when Qu < 1, transmitted 
fiux is possible only through bubbles produced by halos of mass at least ~ 10^^-IO^^Mq. This 
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permits a rough estimate of Agp in this case. Indeed, if C = 0.25 then at z = 6 the number density 
of active bubbles produced by halos with M > W^^Mq is ~ 0.1 per Mpc'^. Each bubble has an 
''max ~ 0.4 Mpc, so if cach bubble produces a low optical depth out to a projected radius of rinax/2, 
then Agp ~ 80 Mpc. This explains the high values of Aqp found in the next section whenever 

Qu < 1. 

3. Results 

To compare our model with the observations, we first note the sizes of the widest dark regions 
found in the quasar spectra. Becker et al. (2001) found a dark region at z = 5.87-6.15, or 16.8 Mpc 

(considering just Lya, not Ly/3), and Djorgovski et al. (2001) found a dark region aX z = 5.25-5.31, 
or 4.7 Mpc. In both cases, the highest transmission spikes within these regions correspond to r ~ 2, 
but measurement noise may well explain these spikes. More quantitatively, within the dark region 
Djorgovski et al. (2001) estimate a mean flux level of 1.6 x 10~^ relative to the quasar continuum, 
and a standard deviation of ~ 0.01 per resolution element (which is ~ 10 pixels). Since there are 
~ 50 resolution elements within the dark interval, both of these numbers imply that the upper 
limit on the flux within a single resolution element corresponds to r ~ 2.5. As our default value we 
adopt this limit of f = 2.5, i.e., we assume that even a small stretch with r < 2.5 can be excluded 
within the dark regions in both of the observed spectra. Figure 1 shows the predicted values of 
Agp and Qu as a function of the cutoff halo circular velocity Vc, for several redshifts between z = 5 
and z = 10. For our standard parameters {Vc = 16.5 km/s at z = 6), we find Agp = 13.3 Mpc and 
Qu = 5.0. At a given redshift, Qu declines more steeply with Vc at high values of Vc, because the 
halo mass function declines steeply at high mass values which correspond to very rare halos. At 
low Vc, Agp depends only weakly on the cutoff since even when Qu > 1 the total cross-section of 
low-T regions in H II bubbles is dominated by bubbles due to relatively massive halos. When the 
minimum Vc is increased, as Qu drops below 1 there is a steep rise in Agp due to the effects of the 
damping wings of the neutral IGM. 

Figure 2 shows the dependence of Agp and Qu on A,- and on f. In every case we set (A) = 2A7-, 
except that (A) is not increased past unity. In general, a lower (A) makes it easier to reach the 
end of overlap, since a greater mass fraction of the gas remains in neutral clumps which do not 
need to be ionized in order to reach Qu = 1- A lower A,- lowers Trcson at a given radius in a 
bubble in proportion to A^, where one factor of A^ is due to the reduced neutral fraction. The 
plot of Agp as a function of f indicates that a more detailed comparison should be possible in the 
future between theory and observation. Indeed, the full probability distribution of r as a function 
of redshift, and the correlation function of r along the line of sight, should be compared between 
detailed observations and detailed numerical simulations with accurate radiative transfer. 

Figure 3 shows the dependence of Agp and Qu on the source parameters Njon and C. Increasing 
A/ion causes a proportional increase in Qu, since a greater number of ionizing photons is produced. 
At the same time, Aqp decreases since the ionizing intensity increases at a given distance from 
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Fig. 1. — The mean free path Agp and the filling factor Qu versus the minimum circular velocity 
Vc for halos which host galaxies. The bottom panel shows the mean free path for crossing through 
a region with a Lya optical depth less than f = 2.5. The mean free path in proper Mpc is shown 
versus the cutoff, minimum halo circular velocity Vc in km/s. The upper panel shows the H II filling 
factor versus Vc- In both panels, the different curves correspond to redshifts 5, 6, 7, 8, 9, and 10, 
as indicated. Vertical dotted lines show the values of Vc that correspond roughly to H2 cooling (2.0 
km/s), atomic cooling (16.5 km/s), and pressure suppression (41 km/s). Other parameter values 
are = 0.4, (A) = 0.8, iVjon = 20, and C = 0.25. 
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Fig. 2. — The mean free path Aqp and the filhng factor Qn. Bottom panels show the mean free path 
in proper Mpc for crossing through a region with r < f. Upper panels show the H II filling factor. 
In the left panels, the gas density A,- is varied, where A = p/p. In every case we set (A) = 2At-, 
except that (A) is not increased past unity. The vertical dotted line shows our standard value of 
At- = 0.4. In the right panels, the maximum optical depth f is varied. The vertical dotted line 
shows our standard value of f = 2.5. Other parameter values are z = 6, Vc = 16.5 km/s, A^ion = 20, 
and C = 0.25. 
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a source. As before, Aqp increases sharply if Nion is lowered to the point where Qu < 1. The 
duty cycle does not affect Qu since we assume that the bubbles of dead sources remain ionized. 
Each active bubble is more highly ionized when <^ is small, but only bubbles with active sources 
can have a low optical depth, so a low ( leads to a slight increase in Agp- The dependence on 
is weak also when parameters are chosen so that Qu < 1. Given this weak dependence, Qu and 
Agp are also expected to be insensitive to clustering of sources. Consider a simple toy model for 
clustering in groups of iV > 1 halos, where we make each source A'^ times brighter but assume 
N times fewer sources. As before, Qu remains unchanged, but each bubble becomes bigger and 
more highly ionized at a given radius. Indeed, the radius at which we find a given Tj-eson goes up 
as VN, so the cross-section goes as oc A^, which cancels out the 1/A decrease in the density 
of sources. Thus, when Qu > 1 clustering has almost no effect on Aqp, but it does decrease Aqp 
somewhat when Qu < 1, since damping wings have a smaller effect on larger bubbles. Note that 
even parameters that do not change \qp significantly can be probed with detailed observations 
since they do change ionized gas fractions within H II bubbles as well as the distribution of bubble 
sizes. 

The discussion so far shows that after considering a wide range of parameter values we expect 
the universe to be well past the end of overlap by redshift six. However, present observations place 
only weak direct constraints on the parameters controlling the abundance and efficiency of sources 
at such high redshifts. Thus, our goal is to use directly the observations of wide regions of Gunn- 
Peterson absorption to infer the state of reionization at z = 6, regardless of the precise parameter 
values. To this end, Figure 4 shows a scatter plot of points derived from the entire parameter 
space of reasonable values. Shown are 5000 points, where each corresponds to a set of parameters 
selected randomly from the ranges z = 5-10, Vc = 1-200 km/s. A,- = 0.2-1, Aion=2-200, and 
^ = 0.01-1. We hold f = 2.5 fixed, and choose parameters from uniform distributions in the log 
(except for z). In the figure we also mark with an x the location corresponding to our standard 
parameter values, and we draw a straight line at Qu = 1- The figure clearly illustrates that over a 
broad range of possible parameter values, there is a very strong correlation between Qu and Acp. 
Indeed, for most parameter values we find that Xqp < 100 Mpc if and only if Qu > 1- Thus, when 
Qu < 1 the incidence of transparent regions should become very low, since for example z = 6-7 
corresponds to a line-of-sight length of only 51 Mpc. Thus, e.g., if Agp = 1000 Mpc then the chance 
of encountering any r < 2.5 spot along a given line of sight is about 1 in 20. Note also that while 
the y-axis in the figure is shown only up to Aqp = 10^ Mpc, most of the points at Qu < 1 actually 
lie at even higher values of Agp- 

The appearance of Figure 4 is determined by the various input parameters whose individual 
effects were shown in the previous figures. If Vc is varied from small to large values, with all other 
parameters fixed, the result is a single curve which goes from the top left to the bottom right corner 
in the figure. Siich curves drawn at various redshifts all lie almost on top of each other at Qu > 1, 
but at Qu < 1 the redshift dependence explains almost all of the scatter in the figure; in fact, if 
we used only z > 6 instead of z > 5, there would be significantly less scatter at Qu < 1, with the 
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Fig. 3. — The mean free path Agp and the filhng factor Qn. Bottom panels show the mean free 
path in proper Mpc for crossing through a region with a Lya optical depth less than f = 2.5. 
Upper panels show the H II filling factor. In the left panels, the source efficiency A^ion is varied. 
The vertical dotted line shows our standard value of A'^ion = 20. In the right panels, the source duty 
cycle C is varied. The vertical dotted line shows our standard value of = 0.25. Other parameter 
values are z = 6, Vc = 16.5 km/s, A,- = 0.4, and (A) = 0.8. 
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Fig. 4. — The mean free path Aqp versus the fihing factor Qn. The vertical axis is the mean free 
path in proper Mpc for crossing through a region with a Lya optical depth less than f = 2.5. The 
horizontal axis shows the H II filling factor. The scattered points show the results from 5000 sets 
of parameters selected randomly from the parameter ranges shown in figures 1, 2, and 3 (except 
that f = 2.5 is held fixed). The large x marks the point corresponding to our standard parameter 
values. The solid vertical line is drawn at Qu = 1. 
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points sticking closer to the vertical line Qu = 1. In order to be conservative we also varied A,- 

over a wide range, and this parameter is responsible for most of the scatter at Qu > 1. As for 
the other parameters, variations in the source efficiency A^ion can be responsible for only a minor 
fraction of the apparent scatter; finally, varying (" has little effect on the figure. 

By any reasonable definition of prc-overlap, bubbles of transmitted flux should be extremely 
rare during pre-overlap. For a given filling factor Qu < 1, the probability that a given point 
is inside at least one bubble is Qu, if the point is chosen uniformly over volume. If we make the 
approximation that different bubbles are independently placed then the probability of a given point 
being in exactly one bubble is Qu exp[— Qn]- Thus, the probability that a given point is in two or 
more overlapping bubbles is 

Pover = Qll [1 - e-«"] . (17) 

Pre-overlap can be defined as the time when Pover ^ 1- For example, Pover = 0.1 when Qu = 0.34 
and Pover = 0.01 when Qu = 0.10. 

4. Discussion 

In the previous section we showed that over a broad range of possible parameters that char- 
acterize the IGM and the population of ionizing sources, there is a strong correlation between the 
reionization state of the IGM as measured by Qu and the typical observed length of Gunn-Peterson 
absorption Agp- We now consider the recent spectral observations in light of these results. Both 
the 16.8 Mpc dark region at 2; = 6 (Becker et al. 2001) and the 4.7 Mpc dark region ai z = 5.3 
(Djorgovski et al. 2001) most likely correspond to Qu > 1- If indeed true, this means that these 
observations correspond to post-overlap just like observations at lower redshift, but their novelty 
comes from the approach toward overlap which implies that the voids, while being mostly ionized, 
still contain a high enough neutral fraction to produce wide regions of high optical depth in spectra. 
The 2; = 5.3 observation cannot correspond to Qu < 1, since there are regions of high transmitted 
flux at still higher redshifts in the same spectrum, while Qu is expected to increase with time (see 
Figure 1), and Qu < 1 should correspond to a spectrum that is almost completely dark except 
for the proximity effect of the quasar (see Figure 4). Note also that we have selected the 4.7 Mpc 
region at z = 5.3 as the widest dark region over the range z = 5.2-5.6, which measures 30 Mpc. 
This selection effect means that the length of 4.7 Mpc must be a few times larger than the effective 
Agp over this redshift range. In contrast, the z = 6 observation may in fact imply the discovery of 
the IGM during overlap, since the dark region is observed right up to the region affected directly 
by the quasar emission. Clearly, however, a dark region of ~ 50 Mpc must be observed before we 
can conclude that we are likely observing the effects of damping wings and of the Qu < 1 era. 
Alternatively, if 15 Mpc dark regions are observed at 2; = 6 along many lines of sight, with the dark 
regions always extending right up to the zone of the proximity effect of the quasar, this also will 
constitute strong evidence that the end of overlap has been observed. 
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We now consider an additional argument which deals with an aspect of reionization that has 
not been included thus far. We examine the effect of the increasing intensity due to multiple 
sources reaching every point within overlapping bubbles. This also has implications for the amount 
of time after the end of overlap during which wide regions of Gunn-Peterson absorption can still 
be observed. We embark on this discussion in order to understand an additional implication of 
the recent observations, which is a result of the significant differences between the spectra along 
the two lines of sight. We argue that the observations are inconsistent with a rapid, homogeneous 
overlap era, and that this presents an important challenge to theoretical models. 

The end of overlap (Qn = 1) niarks the first time when all gas in voids can see at least one 
ionizing soiucc. Thus, as long as Qu < 1, some voids are still completely neutral and they block 
the radiation from some ionizing sources. The last stages of overlap should occur within a small 
fraction of a Hubble time at the overlap redshift, since each source is required to ionize only the 
region out to its nearest neighbors, which are likely to be fairly nearby. In addition, the photons 
which gradually succeed in ionizing the dense surroundings of each galaxy during pre-overlap then 
find it very easy to ionize the low-density voids during overlap. Each time two or more bubbles are 
joined, the gas within the combined volume is exposed to the ionizing radiation from more then one 
source. Therefore, the ionizing intensity inside H II regions begins to rise, allowing those regions 
to expand into high-density gas which had previously recombined fast enough to remain neutral 
when the ionizing intensity had been low. Since each bubble coalescence accelerates the process of 
reionization, the overlap phase has the character of a phase transition and should occur rapidly. 

At the end of overlap, the final shadows cast by neutral gas in voids disappear, and every point 
in the IGM can begin to see all ionizing sources in the universe, except that two barriers must still 
be overcome. One barrier is simply due to the light travel time. At redshift z, the ionizing intensity 
at a given point is determined by taking the number density of sources times the flux of each source, 
and integrating over a spherical volume element out to a radius r. The flux of each source is given 
by its intensity divided by r^. The intensity of each source is proportional to the halo mass, times 
an efficiency factor, divided by the source lifetime. If we use for the purposes of this estimate the 
approximate assumption that all sources have the same lifetime, then the result depends only on 
an integral of halo mass times number density, i.e., on the total collapse fraction -Fcoh which is the 
mass fraction contained in halos which host galaxies. Note that we can set C = 1 since a smaller 
duty cycle increases the intensity of each source but decreases the number density of sources by 
the same factor. For gas of relative density A^, the optical depth which results from integrating 
sources out to radius r is 

-0 78 fl±l\^ f^mh^Y^' fnbff_ ts/C \ [N^F^Y^ [ArY ( r Y^ 
rreson-y.^»(^ 7 ) Y-IA^ ) Vo.0211 5 X lO^ yry* V 20 0.15y' Vo.4; UW ' 

(18) 

The result that the optical depth approaches when r is large arises from the proportionality 
of the total intensity to the radius r. This divergence of intensity with distance is well known from 
Olber's paradox, which applies here since we are considering a small region and ignoring redshift 
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effects. This simple estimate implies that if the light travel time were the only barrier, gas in voids 
would reach r < 2.5 in 13 Myr, and r < 1 in 32 Myr. These numbers can be increased by making 
sources less efficient (i.e., decreasing A'^ion) or by making high-mass halos dominate over low-mass 
ones. As mentioned in §2.1, one natural process which increases Vc is the suppression of gas infall 
into low-mass halos. This suppression may already be important at the end of overlap, and if it 
raises Vc from 16.5 km/s to 41 km/s then Fcoi at 2; = 6 decreases from 0.15 to 0.066, increasing the 
light travel time needed to reach r < 2.5 by a factor of 2.3. If future observations discover that star 
formation at high redshift occurs efficiently only in galaxies with Vc = 100 km/s or greater, then 
the light travel time will be increased by an additional factor of 4.7 beyond the effects of pressure 
suppression. In addition, although we have used the fiducial value = 0.4 in the estimate of 
the optical depth, in this context we are particularly sensitive to voids which have a relatively low 
underdensity and therefore still produce wide regions of Gunn-Peterson flux. If we find a void with 
= 0.8, then in this region the light travel time is increased by a factor of 4. 

The second barrier to increasing the ionizing intensity comes from the shadow due to the small 
volume of gas which remains in neutral clumps even during post-overlap. As in the discussion at 
the beginning of §2.2, we adopt here the model of Miralda-Escude et al. (2000) for the density 
distribution in the IGM and also for the resulting mean free path of ionizing photons. In this 
model, if gas at z = 6 is ionized up to a density of 10 times the cosmic mean then the ionized gas 
has a clumping factor of C = 1.6 times that for gas at the cosmic mean density, and the mean 
free path is 2.9 Mpc. If the gas is instead ionized up to 100 times the cosmic mean density then 
C = 4.5 and the mean free path is 33 Mpc. Thus, the remaining shadow of neutral regions may 
temporarily prevent the intensity from increasing rapidly, at least until most of the gas in filaments 
with moderate overdensities is ionized. A different way to get a similar result is to extrapolate to 
high redshift the absorption due to the Lya forest; Becker et al. (2001) find that this implies an 
overall average optical depth of ~ 2-3 at z = 5.3, and this fairly high value implies a shadow that 
must keep the ionizing intensity from increasing rapidly. 

An important role in delaying the end of the Gunn-Peterson trough may be played by the 
population of photo-evaporating halos. If molecular hydrogen is dissociated, as discussed above, 
then galaxies around the time of overlap form only in halos above ~ 10^ Mq, where atomic cooling 
is efficient. However, most of the gas which lies in halos is inside smaller halos (down to ~ 10^^ Mq), 
where the lack of cooling prevents the formation of galactic disks and stars or mini-quasars. Barkana 
& Loeb (1999) showed that photoionization heating by the cosmic UV background (or shock heating: 
Gen 2001) could then evaporate much of this gas back into the IGM, with the process beginning 
during overlap. They showed that this process affects a broad range of halo masses, with only a 
small gas fraction evaporating out of 10^ Af^ halos, but with halos below ~ 10^ Mq losing their entire 
gas content because of their shallow gravitational potential wells. For overlap at z ~ 6, around 20% 
of the gas in the universe undergoes this photo-evaporation process, and this gas represents most 
of the gas which is already in halos just before overlap. Gas in the smallest halos evaporates very 
quickly into the IGM, but the larger halos which retain some of their gas expand only gradually. 
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and may cast a shadow for some time, preventing ionizing photons from reaching large distances. 
These surviving hales may also contribute to the high column density end of the Lya forest at 
2 ~ 3 (Bond, Szalay, k Silk 1988; Abel k Mo 1998). 

The process of reionization may be very inhomogencous, with overlap occurring at substantially 
different times in different places. Such inhomogeneity can be enhanced due to variations in the 
IGM density, clustering of sources, and spatial or temporal variation in source parameters such 
as the efficiency A^ion- Still, a large inhomogeneity is unlikely if the relatively common low-mass 
halos dominate. Taking the cooling mass as the lower mass limit at 2; = 6, for instance, there are 
3000 halos per Mpc^. The number density is 130 per Mpc^ for = 41 km/s, and 3 per Mpc^ for 
Vc = 100 km/s. Thus, if the efficiency Aion is high only in high-mass halos, then strong Gunn- 
Peterson absorption can last down to lower redshift, because there are fewer stars in galaxies (i.e., 
Fcoi is smaller), because reionization is more inhomogencous, and because the bubble around each 
individual source must grow very large and thus the source must overcome a strong shadow due to 
neutral clumps within the bubble. Whatever the explanation, the two recent observations together 
imply that the increase in the ionizing intensity during post-overlap must have been gradual, or 
that overlap must have occurred inhomogeneously. This follows from the fact that overlap in one 
line of sight occurs at z > 5.9, while a relatively large region of strong absorption is found along a 
second line of sight at z = 5.3, and the age difference between these redshifts is large (140 Myr). 
This significant cosmic variance, which has been observed even with only two lines of sight, further 
justifies the statistical approach that we adopted in §2.3, where we opted to study Agp rather than 
the mean optical depth. 

Note that the the ionizing intensity starts to increase in overlapping bubbles even when Q < 1, 
but this effect should be significant only in the short period when Q approaches unity. Furthermore, 
when Q is significantly less than unity any transmitted flux must be very rare since the damping 
wings block out flux from most bubbles very effectively, almost regardless of the source parameters. 

5. Conclusions 

We have constructed a model of reionization in order to interpret recent observations of strong 
absorption in quasar spectra. By taking the H II bubble produced by each ionizing source, and 
considering the statistics of the ensemble of these bubbles, we have derived for given source param- 
eters the filling factor Qu of H II in the IGM, and the mean free path Agp for observing a region 
with optical depth less than a given f. In particular, we expect the universe to be well past the 
end of overlap by redshift six, with a Agp ~ 10 Mpc. However, regardless of the source parameters, 
we find a strong correlation between the reionization state of the IGM as measured by Qu and 
the typical observed length of Gunn- Peterson absorption Agp. We conclude that the post-overlap 
era is consistent with the observed dark regions of 16.8 Mpc at 2; = 6 (Becker et al. 2001) and 4.7 
Mpc at z = 5.3 (Djorgovski et al. 2001). The z = 5.3 observation cannot correspond to Qu < 1, 
since the universe before the end of overlap should correspond to a spectrum that is almost com- 
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pletely dark except for the proximity effect of the quasar. This strong absorption results from the 
damping wings of neutral IGM, which block out the flux from all except very large, and rare, H II 
bubbles. The z = 6 observation is consistent with the state of the IGM before the end of overlap, 
since the dark region extends right up to the highest available redshift. If this interpretation is the 
correct one then much wider dark regions are expected and should be revealed by spectra of even 
higher-redshift sources. Alternatively, the discovery of the reionization era can be established even 
with quasars at z > 6 if multiple lines of sight all show Gunn-Peterson absorption. 

Combining the two recent observations constrains the evolution of the ionizing intensity during 
and after overlap. The observations imply that along one line of sight, overlap occurs at z > 5.9, 
while a relatively large region of strong absorption is found along a second line of sight at z = 5.3. 
The age difference between these redshifts is 140 Myr, which is much longer than our estimate 
of the light travel time needed to make the voids transparent due to the combined intensity of 
multiple ionizing sources. Thus, the ionizing intensity during post-overlap must have increased 
gradually, with a delay caused by the shadows due to the remaining neutral clumps. The neutral 
gas in halos which photo-evaporated during and after overlap may also have cast a significant 
shadow. If reionization was dominated by relatively high-mass, rare halos, then overlap occurred 
inhomogeneously, and this can also help to explain the observations. 

In this paper we have planted the initial seeds of a very rich phenomenological field. Observers 
and theorists should calculate and compare the probability distribution of r as a function of redshift, 
and the correlation function of r along the line of sight. These statistical quantities will contain 
information on the size distribution of H II bubbles, on the neutral fraction distribution within the 
bubbles, and on the density distribution of neutral IGM outside the bubbles. This emerging field 
will benefit from measurements of increasing resolution and signal to noise, made toward many 
more lines of sight, and extended to higher redshifts. The results in the near future should be a 
definitive detection of the different stages of the reionization era, and the accumulation of a large 
body of knowledge on the properties of ionizing sources and of the IGM during reionization. 

I thank Zoltan Haiman for useful discussions which motivated this work, and George Djorgovski 
for helpful details regarding his observations. I thank CITA for providing research support and a 
stimulating research environment. 
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